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■ Abstract 

> ; 

■ We argue that the knee in the cosmic ray spectrum at energies E ^ 10^^'^ eV is due 
7—1 , to "new physics" , namely to a channel in the high energy ( ^ TeV in the CM) proton 

I interactions hitherto unaccounted for in estimating the energies of the air shower 

' cosmic rays. The new interaction transfers part of the primary particle's energy 

. to modes which do not trigger the experimental arragement (neutrinos, lightest 

I supersymmetric particle, gravitons ) thus underestimating its true energy. We show 

Qh' that this underestimate leads naturally to the observed break (the "knee") in the 

Q . inferred cosmic ray spectrum. The suggestion we advance fits nicely to current 

|H I theoretical extensions of the Standard Model (super symmetry, technicolor, low scale 

^ ' gravity) where new physics at the TeV scale manifests with the distinct signature 

^ ' of missing energy. We present a simple model where the new physics proceeds via 

• i-H , gluon fusion and assuming a single power law for the galactic {E ^ 10^^'^ eV) cosmic 

^ ' ray spectrum, we produce a good fit to the data in the 10^^ — 10^^'^ eV range. Our 

■ proposal should be testable in laboratory experiments (LHC) in the near future 
and, should it be proven correct, it would signal besides the presence of new physics 
in high energy interactions, a drastically different interpretation of the sources and 
acceleration of cosmic rays. 
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The origin of cosmic rays is a subject which, despite the observational and 
theoretical progress made since their discovery has not been settled as yet. 
The reason can be traced to the breadth of their spectrum which extends over 
11 orders of mag nitude to > 10^° eV (see e.g. [1] for a recent review) and 
the fact that by virtue of their diffusion through the galaxy most information 
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concerning their sources is practically lost. Thus, to zeroth order, especially at 
higher energies (> 10^^ eV) at which cosmic ray composition measurements are 
difficult, the sole source of clues about the cosmic ray origin and acceleration 
is their over all spectrum. 

The cosmic ray spectrum consists, roughly speaking, of three distinct sections, 
each of power law form, E'"* in the particle energy but with different values 
for the index 7: In the range 10^ — 10^^'^ eV the index 7 ~ 2.75. Above this 
energy (the "knee"), the spectrum steepens to a power law of 7 ~ 3 which 
extends to ~ 10^^ eV, with some evidence for a further steepening in the 
spectrum indicating a possible cut-off at ~ 10^^'^ eV [2]. This steepening is 
reversed at slightly larger energies (at the "ankle" ) with the spectrum flatten- 
ing to 7 ~ 2 — 2.5 and extending io E 10^°'^ eV, at which point the existing 
statistics are too poor to provide a well defined fiux measurement. 

Considering that cosmic rays propagate in the galaxy by diffusion through the 
tangled interstellar magnetic field, one can argue convincingly that particles 
with gyroradii larger than the galactic scale height (~ 1 kpc) ought to be 
extragalactic. Given that the gyroradius of a proton of energy E[eV) is Rg ~ 
Ikpc Eis/B_Q (where E'lg = £'(eV)/10^^ and B^q is the galactic magnetic 
field in /iG), it is expected that protons of energy E ^ 10^^'^ would escape 
freely from the galaxy. This notion is in agreement with the indication of an 
additional steepening or a potential cut-off in the spectrum at E 10^^ eV, 
with the subsequent fiattening at higher energies being naturally interpreted 
as due to a "harder" extragalactic component. 

This latter extragalactic component, which includes several events above 10^*^ 
eV, has caught recently the attention of the community: It is well known [3] 
that protons of energies ^ 10^^'^ eV suffer catastrophic photopion production 
losses on the Cosmic Microwave Background (CMB). If this extragalactic com- 
ponent permeates uniformly all space, as it was thought to be the case, this 
process should then lead to a cut-off (the so-called GZK cut-off) rather than 
an excess flux above this energy. The potential identiflcation of the source of 
this component with either gamma ray bursts within 100 Mpc [4], a novel, 
neutral hadron immune to the photopion losses [5], or the decay of heavier 
Big-Bang relics [6] has provided the impetus for a recent flare of activity re- 
garding the origin of this speciflc part of the cosmic ray spectrum (see [7] for 
a review). 

However, it is not only this highest energy regime which defles our understand- 
ing of the cosmic ray spectrum origin. The spectrum at energies E ^ lO-*^^"^ 

eV, thought to be galactic (see however [8]) and presumably easier to com- 
prehend, challenges on its own right our understanding of its origin, arguably 
more severely than the extragalactic component, at least from the purely "flt- 
ting the curve" point of view. 
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This last statement, while counter intuitive at first glance, it is nonetheless 
true, as shown by the following simple argument (see also [9]): It is easy to 
obtain a "flattening" of the (any) spectrum by combining two independent 
components, since the harder one will always dominate at sufficiently high en- 
ergies. This appears to be the case with the cosmic ray spectrum at £■ ^ 10^^"^ 
eV. On the other hand, producing a continuous spectrum with a steepening 
break similar to that observed at the cosmic ray spectrum "knee" is much 
harder: It demands the presence of two distinct acceleration mechanisms, one 
of which carries the particles to the "knee" with spectrum oc E~'^-'^^ and a sec- 
ond one which takes practically all the particles that reach the "knee" via the 
first mechanism and only these, to a thousandfold higher energy with spec- 
trum ~ E~^. If this second acceleration mechanism accelerated only a fraction 
of the particles that reach the "knee" (a perfectly "reasonable" assumption for 
most acceleration processes), it would lead to a (not observed) discontinuity 
in the spectrum at this energy. 

To complicate matters further, the most promising acceleration mechanism 
of galactic cosmic rays, namely supernova shocks, can barely produce (even 
theoretically) particles of energies as high as Em — Z ■ 10^^ eV [10] (the 
nuclear charge Z appears because of the dependence of the gyroradius on it) , 
even with the diffusion coefficient at the Bohm value [11]. Energies as high 
as that of the "knee" can be achieved only by assuming that the cosmic ray 
composition at this point consists mainly of Fe nuclei. There exists no known 
(to the authors) mechanism which would carry even a fraction of the (diffusing 
through interstellar space) particles of the "knee" to the energy of the "ankle" , 
in a way that produces the observed spectrum. Furthermore, it is hard to argue 
that the transition between the two power laws at the "knee" could be due 
to some unknown transport effect of the cosmic rays in the galaxy, because it 
appears to be much too sharp for such an interpretation. 

Motivated by the above considerations we are led to propose that the break 
at the "knee" of the cosmic ray spectrum is indicative, not of a distinct accel- 
eration mechanism, but of the emergence of "new physics" in the high energy 
proton interactions, namely of a new channel beyond those considered in the 
models employed to infer the primary particle energy in the air shower arrays. 
We argue that if a fraction of the energy associated with this new channel 
is in a form that docs not trigger these detectors, it will result to an under- 
estimate of the primary particle's energy (a similar suggestion was made in 
[13], based on the interpretation of certain peculiarities in the Extensive Air 
Shower (EAS) data as a sharp change in inelasticity with energy). For a cos- 
mic ray spectrum which is a single power law in energy, this underestimate 
will then manifest as an increase in its slope (a "knee") at the energy at which 
this new channel turns-on, with the spectrum reverting to its original slope 
when it eventually saturates. Furthermore, to account for the break observed 
at the "knee" of the cosmic ray spectrum, this new channel should "turn-on" 
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at an energy ~ TeV at the center of mass, a scale tantalizingly close to that 
at which the emergence of "new physics" is anticipated on the basis of rather 
general considerations [12]. 

As candidates of such new physics we cite supersymmetry (SUSY), a symme- 
try that interrelates bosons and fermions [14] , technicolor, the postulate that 
the Higgs boson is composite made of "techni-fermions" much the way that 
the proton is made of quarks [15], or the postulate that the observed Planck 
scale Mp ~ 10^^ GeV is a 4— dimensional epiphenomenon of the "true" scale 
of gravity which is of order of 1 TeV, "living" in a higher dimensional world of 
D — n + 4 dimensions [16]. According to this last proposal, the gravitational 
interaction would acquire the strength of the rest interactions at the charac- 
teristic scale of ~ 1 TeV, with prolific graviton production (which of course do 
not register in the EAS detectors) in particle interactions above this energy. 
Interestingly, energy which does not register in the EAS detectors (needed to 
account for the break at the "knee" ) is predicted also by the other two propos- 
als: In SUSY this is the hghtest supersymmetric particle (LSP) . In techni-color 
models, pp collisions at energies ^ 1 TeV produce techni-hadrons (hadrons 
carrying techni-color besides the ordinary color); techni-hadron decay leads to 
W's of which a significant fraction decays in neutrinos that also do not trigger 
the EAS detectors. 

One should note from the outset that EAS experiments, by observing the 
results of reprocessing high energy interactions of a largely unknown composi- 
tion through 1000 gm cm~^ of atmosphere, can provide but the most general 
characteristics of these events and as such are not a substitute for the more 
detailed accelerator experiments. However, some of the general properties of 
these interaction can be deduced in these experiments. Such an example is the 
logarithmic increase of the strong interaction cross section with energy, which 
was first inferred from the energy dependence of the height of maximum de- 
velopment of cosmic ray showers. The arguments on the nature of the break 
at the "knee" (i.e. that it is convex rather than concave) are of this general 
character and their interpretation should be viewed in the same spirit. 

To avoid making a specific choice from the list of available alternatives at this 
early stage of our investigation, we model the process simply as the production 
and decay of a system of total invariant mass Mq = 2 TeV and we parameterize 
the entire process by two parameters: the fraction y of the primary particle's 
energy that registers in the cosmic ray detectors and the asymptotic (i.e. 
at energies much higher than the production threshold) ratio a of the cross 
section associated with this new channel to that of the standard interactions. 

On dimensional grounds, the cross section of the new channel is assumed to 
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be of the form 

= ^g{r) (1) 

where B is a dimensionless constant (related to a), s = 2mpE and g{r) is a 
function of the dimensionless ratio r = M^/s. At high energies we expect the 
pp interactions to be dominated by gluon scattering and accordingly 

1 

g{T) = J f{x)f{T/x) dx (2) 



where f{x) is the gluon distribution within the proton, which we parameterize 
as 

/(a;) = i(iV+l)^i^^^ , with iV = 6 (3) 



For the conventional pp interactions the cross section rises slowly with energy 
and for our purposes we consider it to be a constant (7o{E) ~ 80 mbarn. At en- 
ergies above the new physics threshold the cosmic ray interactions will proceed 
either through the standard channels with probability Po{E) — ao{E) /[ao{E)-\- 
(Tn{E)] or through the new channel with probability Pn{E) = an{E)/[(To{E) + 
an{E)]; given that rgi^r) 1/Cq = constant for r ^ 0, setting B = MqCqUoO. 
we get an{E) / ao{E) a for E ^ M'^/2mp as desired. 

Whenever high energy cosmic ray particles interact through the new channel, 
events of total energy E' will register at the detector an energy E = yE' {y < 
1). Therefore, if the cosmic ray intensity is I{E), for particles interacting 
through this new channel, the inferred intensity will be of the form 

J I{E')Pn{E')5{E - yE') dE'^^I j P. j (4) 

while for events interacting through the conventional channel the resulting 
intensity will have a similar form but with y — 1 and Po{E) in place of 
Pn{E/y). 

Assuming the incident galactic cosmic ray spectrum to be of the form Ijg{E) = 
E'^^exp^—E / Eq) with 7 ~ 2.75 and Eq ~ 10^^'^ eV, a value consistent with 
the earlier argument on the cosmic ray gyroradii at E ~ Eq, the observed 
cosmic ray flux at an energy E will be 

^ h + Cargir)'^ l/[C,{Ty)g{Ty)]+a\ 
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The first term in the sqiiarc brackets in Eq. (5) represents the contribution 
to the spectrum from interactions through the conventional channels while 
the second that due to the new one. The presence of the exponential cut-off 
broadens and deepens the effects of the presence of the new channel. Their 
combined effect is necessary for a good fit to the data. 
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Fig. 1. The cosmic ray spectrum f{E) multipHed by E'^ for E > 10^^ cV. Long 
and short dashed lines are respectively the incident galactic and extragaiactic 
components. Solid line is the combined spectrum expected to be measured for 
y = 1/2, a = 2. +'s are the Tibet data, diamonds the Casa Bianca data, squares 
the Fly's Eye data and x 's the AGASA data. 

In addition to the galactic cosmic ray component, whose contribution is ex- 
pected to be unimportant beyond Eq ~ lO^^eV, there exists also an extra- 
galactic component, whose reprocessing in the atmosphere should also re- 
sult in a modification of its spectrum according to the prescription of Eq. 
(5). The precise form of this component is of course unknown since it is 
dominated at lower energies by the galactic component. Following [1] we 
assume its spectrum to be of the form Ieg{E) oc E^'^exp{—E/Ei) with 
5 ~ 2.2 - 2.6 and ~ 10^° - lO^"-^ eV. In Figure 1 we present the en- 
tire (galactic + extragaiactic) cosmic ray spectrum from 10^^ — 10^^ eV, with 
7 = 2.75, g = 2.2,^0 = 2 lO^^eV, = lO^O-^eV by applying the effects 
of the new postulated channel in the interaction to both components with 
y — 1/2, a — 2. We also plot the relevant data from two different experiments 
in each of the lO^^-lO^*' eV (Tibet [17], CASA-BLANCA [18]) and lO^^-lO^o 
eV (Fly's Eye Stereo [2], AGASA [19]) energy ranges. We expect that these 
should bracket the true values of the corresponding parameters and should 
serve as a gauge of the systematic errors involved in computing the cosmic ray 
spectra in each range. 
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While it is very difficult to draw immediate conclusions favoring specific mod- 
els from the existing data our general considerations appear to be on the 
correct footing: Our proposal produces a break in a single power law cosmic 
ray spectrum of the correct character (an increase in the slope) over a very 
limited range in energy (half a decade), as demanded by the data. As noted 
earlier, both these features are very hard to produce with more conventional 
approaches. Our proposal is further supported by the more detailed analysis of 
EAS data: the apparent very sharp change in cosmic ray composition to almost 
exclusively Fe, inferred from the abrupt change in the depth of the maximum 
in the shower development around 10^^ eV (fig.5 of [18]), is quahtatively of the 
form expected by a sharp increase in the interaction cross section, such as we 
propose, and an ensuing dispersion of the available energy to a large number 
of secondary particles. Such an abrupt change in composition is very hard to 
understand (let alone model) within the realm of conventional models. 

Considering the fit of our figure at energies higher than the "knee" , one notes 
that the assumption of a cut-off in the galactic component at ~ 19^^'^ eV is 
necessary for a good fit, but it is also reasonable, supported both theoretically 
(the gyroradii arguments above) and experimentally (as discussed in [1]) by 
the observed increase in the cosmic ray anisotropy at this energy [19]. Given 
the simplicity of our assumptions we think that this fit is particularly good. 
One could think of several ways of improving it, if necessary, at the expense of 
introducing additional parameters in modeling the physics of the high energy 
interactions. We think that the present quality of the data does not warrant 
at this point such an extension of the assumptions used. The presence of the 
extragalactic component does affect the values of the fitting parameters since 
this component does make some contributions at lower energies too (the values 
of q and Eq are thus closely related). However, this seems to affect little the 
values of y and a which determine primarily the fit of the spectrum at the 
"knee". 

Where all these leave us? Our interpretation carries with it a number of con- 
sequences: (a) To start with, it implies the presence of "new" structure in the 
high energy physics interactions at energies consistent with those suspected on 
the basis of generic theoretical considerations. The new physics is shghtly be- 
yond the reach of the Fermilab Tevatron, but it will be preeminently present 
at LHC Their "benchmark" signature, based on our interpretation of the 
"knee", should be a new channel with strong interaction cross section and 
missing energy. We have computed the cross-section of the new channel at 
Tevatron energies using the formalism presented above and it was found to 
be a small fraction of the total cross section. At this point we would like to 
refrain from more specific predictions given the simplicity of our model and 
the quality of the cosmic ray data; however, the above requirements are quite 
strong and may already be in disagreement with some of the above proposed 
alternatives, (b) On the cosmic ray physics side, the demand for a break at the 
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"knee" much sharper than allowed by cosmic ray transport coupled with our 
proposal and fits, lead to the radical suggestion that the cosmic ray sources 
must, by and large, produce spectra extending to the "ankle" (rather than the 
"knee"). This then leads to the unsettling conclusion that supernovae should 
not be the dominant contributor to the cosmic ray spectrum at least at ener- 
gies greater than ~ 10^^ — 10^^ eV. Leaving for the moment the issue of cosmic 
ray composition aside, our interpretation of the cosmic ray spectrum at the 
"knee", implies the presence of an alternate source of cosmic rays capable of 
accelerating particles to the energies of the "ankle". It is interesting to note 
that independent considerations recently pointed to similar conclusions [21]. 
Hints to the nature of these sources may in fact be provided by the observed 
anisotropy at ~ 10^^ eV toward the galactic center [20]. We plan to revisit 
both these issues in a future publication. 

While this paper was being written, the potential effects of physics beyond 
the Standard Model were announced (deviation of the muon g — 2 value from 
that of the standard model [22]). This effect was interpreted as requiring the 
presence of new physics at ~ TeV scales, similar to that involved in our con- 
siderations. 
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